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Summary 

The  above-named  project  was  sponsored  by  the  Naval  Research 
Laboratory  for  a  period  of  two  years,  beginning  on  March  1,  1981,  and 
ending  on  March  1,  1983,  and  was  subsequently  given  a  without-cost 
extension  to  July  31,  1983.  The  greater  part  of  the  research  was 
conducted  by  Professor  Yih-Hsing  Pao  and  Dr.  George  C.C.  Ku  of  Cornell 
University. 

Research  results  are  shewn  in  the  three  technical  reports  included 
herein,  the  titles  of  which  are: 

1.  Acoustic  Radiations  of  a  Summerged  Cylindrical  Shell 

2.  Acoustic  Radiation  from  a  Submerged  Finite  Circular  Cylindrical 
Shell  Driven  by  a  Harmonic  Force 

3.  Attenuation  Constants  and  Phase  Velocities  of  Eighteen 
Polymeric  Materials 

Reference  should  be  made  to  each  abstract,  which  summarizes  the 
results  of  each  investigation. 
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The  project  entitled  "Radiation  and  Diffraction  of  Underwater  Acoustic 
Waves"  was  proposed  by  Professor  Yih-Hsing  Pao  of  Cornell  University  in 
November  of  1980  after  consultation  with  Dr.  Lawrence  Flax  of  the  Naval 
Research  Laboratory.  Through  the  Office  of  Naval  Research,  the  project 
was  funded,  beginning  on  March  1,  1981,  for  a  two-year  period.  It  was 
extended  (without  addition  of  funds)  for  seven  months,  ending  on  July 
31,  1983. 

The  Statement  of  Work  of  the  original  proposal  was  revised  in  the 
contract,  with  the  following  two  topics  receiving  focus: 

(1)  Radiation  of  Transient  Pulse  from  an  Elastic  Body 

(2)  Generalized  Ray  and  Geometric  Theory  of  Diffraction. 

To  maintain  close  collaboration  with  Naval  personnel.  Dr.  George  C.C.  Ku 
of  Cornell  University  was  assigned  to  work  at  the  Physical  Acoustics 
Branch  of  the  Naval  Research  Laboratory  under  the  supervision  of  Dr. 

Flax.  Dr.  Flax  left  NRL  in  July  1981,  and  the  research  work  has  changed 
to  concentrate  on  the  analysis  of  radiation  of  a  submerged  circular  cylin¬ 
drical  shell,  under  the  supervision  of  Dr.  L.R.  Dragnette  who  also  directed 
the  experimental  investigation. 

The  analysis  was  based  on  a  mathematical  theory  closely  representing 
the  experimental  model,  and  the  objective  was  to  determine  the  radiation 
pressure  of  the  acoustic  wave  field  near  the  shell  (one  to  ten  radius 
from  the  central  axis  of  the  shell) ,  under  the  excitation  of  a  point 
force  with  a  harmonic  time  function.  Analysis  of  the  motion  of  the  shell 
was  based  on  Donnell's  equation  for  circular  cylindrical  shells,  and 
that  of  the  acoustic  radiation  was  based  on  the  scalar  wave  equation. 


(iii) 


The  shell  vibration  is  coupled  to  the  fluid  particle  motion  through  the 
reactive  fluid  pressure  on  the  surface  of  the  shell. 

Both  shells  of  infinite  length  (neglecting  the  reflections  of  waves 
from  both  ends  of  the  shell)  and  finite  length  have  been  investigated. 

The  results  for  an  infinite  shell  are  presented  in  Technical  Report  No. 

1,  which  is  included  with  this  report.  The  analysis  follows  the  previous 
work  by  W.  Voger  and  D.  Feit,  "Response  of  Point-Excited  Infinitely  Long 
Cylindrical  Shell  Immersed  in  Acoustic  Medium  (NSRDC  Report  80-061,  1980) . 

Der tailed  angular  variations  of  the  radiation  pressure  near  the  shell 
are  calculated  at  several  axial  locations.  These  results  were  compared 
with  experimental  measurements  which  are  not  included  in  this  report. 

Analysis  of  the  radiation  for  a  finite  cylindrical  shell  is  discussed 
in  Technical  Report  No.  2.  The  problem  is  far  more  complicated  than 
that  for  an  infinite  shell,  and  numerical  calculations  have  not  yet  been 
completed. 

During  the  second  year,  in  collaboration  with  other  personnel  of 
the  Physical  Acoustics  Branch  of  NRL,  George  Ku  also  worked  in  another  research 
project  on  the  deformation  of  optical  fibers  under  the  radiation  pressure. 

One  phase  of  the  research,  the  measurement  of  phase  velocities  and  atten¬ 
uation  constants  of  eighteen  polymeric  materials  was  conducted  at  the 
Ultrasonics  Laboratory  of  Cornell  University  by  Mr.  C.P.  Chen.  The  results 
are  compiled  in  Technical  Report  No.  3,  prepared  by  Mr.  Chen. 

Each  technical  report  contains  an  abstract  which  summarizes  the 
results  of  each  investigation. 


Technical  Report  No.  1* 


ACOUSTIC  RADIATIONS  OF  A  SUBMERGED  CYLINDRICAL  SHWT.T. 


George  C.  C.  Ku  and  Yih-Hsing  Pao 
Department  of  Theoretical  and  Applied  Mechanics 
Cornell  University,  Ithaca,  New  York  14853 


Abstract 


The  acoustic  radiation  of  an  infinitely  long  cylindrical  shell 
submerged  in  water  is  analyzed  on  the  basis  of  the  Donnell's  equation 
for  thin  shell  coupled  with  the  wave  equation  for  the  fluid  pressure 
field.  The  shell  is  excited  into  steady  vibration  by  a  concentrated 
harmonic  force  applied  normal  to  the  shell  surface.  Dispersion  of  waves 
in  the  shell  is  analyzed  and  confirmed  with  the  Pochhammer 1 s  theory  for 
thick-walled  circular  cylinders.  The  comparison  establishes  the  upper 
limit  of  frequencies  for  which  the  analysis  is  applicable.  Radiated 
pressure  fields  around  the  shell  are  calculated  for  a  few  receivers  at 
various  locations  in  the  fluid. 
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INTRODUCTION 


The  acoustic  radiation  of  a  submerged  infinitely  long  cylindrical 
shell  has  been  investigated  by  many  researchers  [1-9].  The  practical 
importance  of  these  investigations  for  the  design  of  underwater  vehicles 
is  well  documented,  but  the  approach  has  been  approximate  in  many  cases. 
One  approach  is  to  apply  Green's  function  of  the  wave  equation  for 
acoustic  pressure  and  an  assumed  vibrational  mode  of  the  submerged 
object  [2-4].  Another  approach  is  to  couple  the  wave  equation  for  an 
infinite  acoustic  medium  to  the  thin  shell  theory  through  the  acoustic 
pressure  on  the  shell  surface  [5-9].  The  mathematical  complexity  of  the 
second  approach  arises  from  the  fluid  reaction,  which  affects  the 
vibration  of  the  shell,  being  itself  determined  by  the  amplitude  of  the 
modes  of  vibration.  A  general  treatment  of  the  interaction  of  acoustic 
waves  with  a  structure  is  given  in  Ref.  10. 

The  problem  treated  in  this  report  is  an  infinitely  long 
cylindrical  shell,  freely  suspended  in  the  fluid  medium  and  vacuous 
inside.  The  shell  is  subjected  to  a  harmonic  point  force  on  the 
cylindrical  surface.  The  wave  propagation  in  the  surrounding  fluid 
medium  is  governed  by  a  simplified  wave  equation.  The  material  of  the 
shell  is  assumed  to  be  isotropic,  homogeneous ,  and  linearly  elastic. 

The  shell  thickness  and  deflections  are  sufficiently  small  that 
Donnell's  equation  for  thin  shells  can  be  applied.  The  dynamic  effect 
of  the  hydrostatic  pressure  is  neglected.  Through  the  condition  of 
continuity  for  the  radial  velocity  of  the  vibrating  shell  and  the  radial 
particle  velocity  of  the  fluid  medium,  the  acoustic  pressure  can  be 
treated  as  a  reactive  force  acting  on  the  submerged  cylindrical  shell. 


The  deflection  of  the  shell  and  the  radiating  acoustic  pressure  are 
described  in  the  cylindrical  coordinate  (r,6  ,z).  Due  to  azimuthal 
symmetry  of  the  loading  and  the  shell  geometry,  the  displacements  are 
conveniently  decomposed  into  the  Fourier  series  in  the  angular 
coordinates  6  .  The  axial  dependence  (in  the  z  coordinate)  is 
eliminated  by  using  the  integral  Fourier  transformation.  Through  the 
condition  of  continuity,  the  acoustically  radiated  pressure  is  derived 
in  terms  of  the  radial  deflection  of  the  vibrating  thin  shell.  To 
calculate  the  radiated  pressure  field,  a  special  numerical  integration 
scheme  is  developed  to  evaluate  the  inverse  Fourier  trams form. 

In  Donnell's  theory,  the  extensional,  torsional,  and  radial 
vibrational  modes  of  a  cylindrical  thin  shell  are  all  coupled  in  a 
system  of  three  equations.  To  show  the  coupling  effect  of  these  modes, 
the  dispersion  curves  for  each  circumferential  mode  are  calculated. 

Three  brandies  of  the  dispersion  curves  are  formed  and  each  one  can  be 
identified  as  one  vibrational  mode  with  its  distinct  group  and  phase 
velocities.  By  comparing  these  dispersion  curves  with  those  derived 
from  the  exact  theory  of  elasticity  for  thick-walled  cylinders  [ll,  12], 
the  validity  of  Donnell's  equations  is  established.  The  agreement  for 
the  fundamental  modes  is  excellent,  but  the  situation  is  somewhat 
different  for  the  higher  modes.  Furthermore,  the  dispersion  curves 
based  on  Donnell's  equation  do  not  exhibit  the  overtones  above  the  first 
three  modes.  This  indicates  that  Donnell's  equations  of  motion  are  not 
adequate  for  evaluating  the  response  of  a  cylindrical  shell  at  high 
frequencies. 


Numerical  examples  are  shown  for  the  radiated  acoustic  pressure 


around  a  circumference  of  radius  equal  to  ten  times  the  mean  shell 
radius.  The  perturbed  pressure  field  (the  perturbation  over  the  ambient 
pressure)  decreases  rapidly  from  the  point  of  application  of  the 
harmonic  force  in  both  circumferential  and  axial  directions.  The 
pressure  varies  in  the  forcing  frequency,  and  it  is  higher  at  100  KHz 
than  that  at  10  KHz  because  the  radial  modes  of  the  shell  are  strongly 
excited  at  high  frequencies. 


2.  EQUATIONS  OF  MOTION  FOR  A  SUBMERGED  CYLINDRICAL  SHELL 

Consider  an  infinitely  long  circular  cylindrical  shell  submerged  in 
a  fluid  medium  as  shown  in  Fig.  1.  The  center  line  of  the  cylindrical 
shell  coincides  with  the  z  axis  of  the  cylindrical  coordinate  (r,  0, 
z).  An  oscillatory  point  force  of  circular  frequency  u>  is  driven  in 
the  radial  direction  against  the  cylindrical  shell.  It  generates  a 
steady-state  motion  in  the  shell  which  radiates  acoustic  energy  into  the 
surrounding  fluid  medium  at  the  same  frequency. 

For  a  thin  cylindrical  shell  of  mean  radius  a,  the  equations  of 
motion  (Donnell's  equations)  for  the  three  components  of  a  displacement 
vector  u(  0 ,z,t)  are 
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where  uf  v  and  w  are  the  longitudinal,  tangential,  and  radial 

displacements  of  the  cyl  indrical  shell  [13,  14].  The  parameter  E  is 

Young's  modulus,  v  the  Poisson  ratio,  Ps  the  mass  density  of  the 
2  2 

shell  and  C  “  E/p  (1-  v  )  where  C  is  the  extensional  wave  speed  in 
P  5  p 

a  thin  plate.  The  last  term  on  the  right-hand  side  of  Eq.  (lc) 
represents  the  harmonic  point  force  placed  at  z*0,  and  5*0  plus  the 


reacting  acoustic  pressure  p(r,8,2,t)  *  pQ(r  ,e  ,z)exp(-i<i)t)  at  the  interface 
between  the  shell  and  the  fluid  medium.  Due  to  the  symmetry  about  a 
diametral  plane,  the  displacements  can  be  expressed  in  the  following 


forms: 

00 

u(0,z,t)  »  \  U  (z)cosn6e  IWt 
n-0 

00 

v(0,z,t)  »  l  V  (z)sinn0e  ltot 
n-0  n 

00 

w(0,z,t)  ■  V  W  (z)cosn0e  IWt 
n-0 

The  index  n  is  called  the  circumferential  modal  number  which  indicates 
a  mode  of  2n  nodal  points  along  the  circumference.  Similarly,  the  delta 
function  in  0  is  expanded  into  a  cosine  series 

5<e)  ■hW05"9 
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2  ,  n  >  1 


Substituting  Bqs.  (2)  -  (4)  into  Eq.  (1)  and  applying  the 
orthogonality  conditions  for  coan9  and  sinn6  ,  we  obtain  the 
following  ordinary  differential  equations  for  ^n(z) ,  Vn(z)  and  Wn(z): 
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These  equations  are  solved  by  applying  the  Fourier  transform,  which 
is  defined  as 


V» 


e  ikzU  (z)dz  ,  etc. 
n 


The  inverse  transform  of  Un(k)  is 


1 

2n 


elkzUn(k)dk 


where  k  is  physically  the  wave  number  along  the  z-axis.  The 
transformed  equations  are: 
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where  3  •  ua/C  is  a  dimensionless  frequency. 
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3.  EQUATION  FOR  ACOUSTIC  WAVES  IN  A  FLUID  MEDIUM 

The  perturbing  pressure  field  P(r,  9  ,z,t)  in  an  ideal  fluid 
medium  caused  by  the  vibrating  shell  is  governed  by  the  wave  equation 
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“F  P  ( f ,  9 , 2 ,  t ) 


where  V2  is  the  Laplacian  operator  and  c  is  the  wave  speed  in  the 
fluid.  For  a  harmonic  time  dependent  wave  motion,  Sq.  (9)  reduces  to 
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where  kQ*  m/c  is  the  wave  number  in  the  infinitely  extended  fluid 
medium.  Applying  again  the  Fourier  transform,  Eq.  (6),  to  the  z- 
var table,  we  obtain 
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This  homogeneous  differential  equation  has  a  solution 


P _( r , Q , k)  «  l  A  H  (<r)cosn9 
u  n  n 


n»0 


where  <  ■  (k2  -  k2}1^2.  The  Hankel  function  of  the  first  kind, 

o 

B  (  <r) ,  is  used  to  represent  outgoing  waves  in  the  infinitely 
n 

extended  fluid  medium.  The  conventional  superscript  (1)  for  Hn  is 
omitted  throughout  this  report. 

The  unknown  coefficients  A  of  Eq.  (12)  can  be  determined  from  the 

n 

continuity  conditions  at  the  shell-fluid  interface.  For  an  ideal  fluid 
the  normal  velocity  must  be  continuous. 
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The  fluid  particle  velocity  vf (a>  q,  2,t)  ia  related  to  the  perturbed 
pressure  P(a,  9  ,z,t)  by  the  momentum  relation 


3vf 

VP(r,8,z,t)  «  -p-j~  (r«9.z»t) 

Under  the  assumption  of  infinitesimal  displacement  of  the  vibrating 
shell,  the  surface  normal  p  of  the  cylindrical  shell  still  orients  in 
the  radial  direction,  and  the  radial  component  of  Eq.  (14)  is 
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Substituting  Eqs.  (2)  and  (12)  into  (IS),  the  coefficient  Ar  is 
obtained  as 

pw2W  (k) 
a  _  _  n 


<Hn(<a) 


Consequently,  the  Pourier  transformed  pressure  field,  Eq.  (12) ,  may  be 
written  as 


PQ(r >9,k)  *  l  Znr (k)Wn (k)cosn0 


where 
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nr 
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On  the  surface  of  the  cylindrical  shell,  the  reacting  acoustic 
pressure  is  written  as 
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4.  RESPONSE  OF  THE  CYLINDRICAL  SHELL 

Replacing  PQ(a,  9,k)  in  Eq.  (8c)  by  that  given  in  Eq.  (19),  we 
arrive  at  the  following  relation: 
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Froa  Bqs.  (8a) ,  (8b)  and  (21) ,  the  radial  displacement  W  (k)  is 

n 

determined  by  using  Cramer's  rule, 

5  00 


where 
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A  -  [ft2-  n2-k2a2]  [ft2-  ^  k2a2  -  n2]  -  k2a2n2 


Note  that  Eq.  (22)  describes  the  relation  between  the  Fourier 
transformed  radial  displacement  W  (k)  of  the  cylindrical  shell  for 

£ 

each  circumferential  mode  n  and  the  harmonic  driving  force  Fexp  (-i'M 

y,  t)  6  { 9  )  6  (z).  Therefore,  Z  (k)  +  Z  (k)  +  Z  (k)  is  the  total 

np  nm  na 

impedance  submerged  cylindrical  shell  per  unit  circuruf erence.  Of 


the  three  tarns,  Z  (k)  is  due  to  the  reacting  fluid  medium,  z  (k) 

na  np 

is  due  to  the  flexural  vibration  of  the  shell,  and  Z  (k)  is  due  to 

nm 

the  extensional  vibration. 

To  determine  the  radial  displacement  f?n(z)  from  (k)  of  Sq. 
(22) ,  the  inverse  integral  Fourier  transform  (7)  is  used.  From  Sq. 

(2c) ,  we  obtain  finally  the  radial  velocity  of  the  submerged  cylindrical 
shell 
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where 


Z(k)  ‘  Znn(k)  +  Z„Jk>  +  Z  <k>  * 
n  np  nm  na 

Similarly,  the  axial  and  tangental  velocities  5  and  v  are  obtained: 
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u  A  2  Z 

Zv(k)  -  a  tk2a2-J  -lyi  +  (-k2a2-  iyi  n2+  a2)] 


Note  that  the  integrations  in  Equations  (24) ,  (26)  and  (27) 
contribute  to  the  velocities  of  the  submerged  thin  shell  for  each 
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particular  circumferential  mode  n,  and  the  total  displacement  is 
obtained  by  summing  up  all  circumferential  mod es.  Por  the  case  of  n-0 , 
8qs.  (24)  and  (26)  show  that  the  radial  and  longitudinal  displacements 
are  independent  of  9.  On  the  other  hand,  Bq.  (27)  indicates  that  the 
tangential  displacement  vanishes.  This  may  be  interpreted  by  the 
decomposition  of  a  point  force  into  an  infinite  series  in  Bq.  (3) .  When 
n»0,  the  point  driving  force  degenerates  into  a  ring  force  at  Z «o  of 
amplitude  P/2ira.  Furthermore,  in  Bqs.  (24),  (26)  and  (27),  if  the 
wave  number  k  in  the  inverse  Fourier  transform  degenerates  to  zero, 
the  radial  and  tangential  displacements  will  be  independent  of  z  and 
the  longitudinal  displacement  vanishes  identically.  This  corresponds  to 
the  case  for  a  cylindrical  shell  with  forces  distributed  along  a  line  in 
the  axial  direction,  normal  to  the  surface  of  the  shell. 


5.  PRESSURE  FIELD  OF  RADIATION 

Th*  radial  notion,  w(  9,z,t)  in  Bq.  (24)  was  derived  in  a  recent 
report  by  Voger  and  Fait  [ is] .  They  also  calculated  the  radial 
displacenents  for  various  orders  of  n.  In  this  report,  we  evaluate  in 
detail  the  radiation  pressure  field  in  the  fluid  medium. 

Substituting  Bq.  (22)  into  (17) ,  the  transformed  acoustic  field  of 
radiation  is  obtained 

—  00  sn^ 

po(r,0,k>  -  l  ;-z-  (t)-z  (k)cosn9  • 
n*0  n 

Thm  pressure  field  is 

P(r.8,2,t)  l  e^cnSe'1"'  f  Znr(k)  |20dk  . 

n»0  *o  n 

Note  that,  in  the  above  equation,  if  the  wave  number  k  <  k  ,  the 

o 

pressure  field  P(r,9,z,t)  is  a  progressive  wave  inasmuch  as  Z  (k)  in 

nr 

Bq.  (18)  acts  as  both  resistive  and  reactive  impedance.  On  the  other 
hand,  when  k  >  kQ,  the  argument  of  the  Hankel  function,  <r  ■  (kg  - 
k  )^  r  becomes  purely  imaginary.  If  the  Sankel  function  of  the  first 
kind  with  the  imaginary  argument  is  replaced  by  the  modified  Bessel 
function  of  the  second  kind, 

Hn(ix)  -  (2/7r)i’(n+1^n{x)  > 


2nr(k)  in  Bq.  (31)  is  given  by 
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ipckoK  [(k2-k£)*r] 

Z(k)  -  - 2 - 7— 

nr  (k2-k2)±K;[(k2-k2)*a] 


It  la  a  pure  reactive  impedance.  Thus  no  acoustic  energy  is  radiated 
froa  the  cylindrical  shell.  Figure  2  shows  the  relation  between  the 
acoustic  wave  number  kQ  in  the  fluid  medium  and  its  projection  along  the 
x  axis  of  the  shell.  It  indicates  that  the  radiation  is  possible  only 
when  k  <  kQ. 


6.  DISPERSION  RELATIONS 

Since  the  displacements  u,v,w  and  the  pressure  field  P  are  all 
related  to  a  coenon  factor  Zn(k)  ■  Znp(k)  *  2^(10  +  zna^'  is 
important  to  examine  the  characteristics  of  this  impedance  for  a 
submerged  cylindrical  shell.  We  consider  first  a  cylindrical  shell  in  a 
vacuum  driven  by  a  point  force. 

Without  the  surrounding  fluid  medium,  the  radial  displacement 
Wn(k)  in  Bq.  (22)  is  reduced  to: 

_  £  F 

Vk>  •  -'IWialz"  W-fi  mTkTT 

np  nm 

The  zeros  of  the  impedance  function  are  given  by  the  roots  of  the 
following  equation: 


2S(k-n>  ■  Zno(k'n)  *  Znr,<k’n>  *  0  • 


np  nm 

When  n-0.  The  above  is  simplified  to: 


z°(k,n) 


r  h2  ,,  \  i* ,  ,  n2,  k2a2v2 

Qa  [  12a2  (ka)  '  Q  £22  -  k2a2 


It  is  a  second  order  equation  in  G  except  when  G  *  oja/Cp  *  0 
There  are  two  roots  for  each  value  of  ka,  and  the  locus  of  the 
roots  forms  two  branches  of  dispersion  curves  as  shown  in  Fig.  3a.  As 
ka  approaches  infinity,  these  two  branches  of  dispersion  curves 
approach  the  following  asymptotes: 


G/ka  -*•  kh//T? 

2-  »  i  + 

.  1 -k2h2/1 2 


as  ka  -*■ 


On  the  other  hand,  when  ka  approaches  zero,  Eq.  (36)  is  reduced  to: 

•  ^{1  -a2+^3'° 

The  branches  of  dispersion  curves  start  at  the  following  values  along 
the  ft -axis: 


‘(38) 


ft  -*•  /I  -  v*  ka  £  0 
a  -  1 


as  ka  +  0 


(39a) 

(39b) 


By  substituting  the  value  of  ft  in  Sq.  (39a)  into  the  homogeneous 
equation  (8) ,  the  eigenfunctions  can  be  easily  identified  as: 


“0<k> 


?o(k)  »  0  as  ka  0 


This  implies  that  the  dispersion  curve  started  at  ft  >0  corresponds  to 
the  extensional  mode  of  a  cylindrical  thin  shell.  Similarly,  if  the 
value  of  ft  >1,  which  corresponds  to  the  ring  frequency  f  ■  0^/(2  ira) 
of  the  cylindrical  shell,  in  Eq.  (39b)  is  used,  the  eigenfunctions  are: 


,  ?Q(k)  *  0  as  ka  -*•  0 


W0<k) 

3 

1 

U0(k) 

vka 

Therefore  the  dispersion  curve  starting  at  ft  *1  corresponds  to  the 
radial  mode  of  the  cylindrical  thin  shell. 

Prom  these  dispersion  curves,  the  phase  velocity  and  the  group 
velocity  for  each  vibrational  mode  are  calculated  as  shown  in  Pig.  3b, 
3c.  Both  the  phase  and  group  velocities  of  the  radial  mode  approach 
Cp,  the  extensional  wave  speed  in  a  thin  plate,  asymptotically.  On  the 


other  hand,  tha  phase  and  group  velocities  of  the  extensional 
vibrational  node  approach  zero  when  lea  is  less  than  10  and  then 
increase  gradually  as  ka  increases. 

For  n*l  in  Eq.  (35) ,  the  dispersion  relation  is: 


Z?(k,n)  »  (a2k2+l  )2  +  1  -  n2]{Cn2-  -  k2a2][a2-  IfW-l] 

12a2  6 

-  k2a2-^)2}+[k2a2(l+vQ2-k2a2v)  +o2-  ^  -  k2a2]  *  0  . 


2 

It  is  a  third  order  equation  in  a  ,  and  there  are  three  distinct 
branches  of  dispersion  curves  for  n-1  node  as  shown  in  Fig.  4a.  When 
ka  approaches  infinity,  Eq.  (42)  becones: 


ZS(M) 


ka^» 


-k*a6  - 


k2h2  .  a 
~TZ~  *  Fa2" 


ft 


3* 


-  1 


ffflsL 

'JipF 


1-vT 

-~r 


This  yields  the  following  asymptotes  for  the  dispersion  curves: 

n/ka  ■*.  kh/vT? 

n/ka  -*■  1  as  ka  » 

n/ka  -  [(1  -  v)/2]* 


Am  ka  approaches  zero,  Eq.  (42)  degenerates  to: 


Z?(M) 


ka-0 


O,  r  - 

=  [3“ 


l]s 


1}  -  0 


This  implies  that  the  three  branches  of  dispersion  curves  should  start 


at: 


ft  *  0 

(46a) 

ft  =  [(1  -  v)/2]*  as  ka  -*■  0 

(46b) 

ft  =  /? 

(46c) 

Again,  the  phase  velocities  and  group  velocities  can  be  calculated  from 
these  three  dispersion  curves  as  shown  in  Fig.  4b  and  4c. 

In  a  case  analogous  to  that  for  n*0 ,  the  phase  and  group 
velocities  of  the  radial  node  (top  branch)  in  Fig.  4a  approach  the 
extensional  wave  speed  of  a  thin  plate.  For  the  extensions!  node 
(bottom  branch)  in  Fig.  4a,  the  phase  and  group  velocities  approach  zero 
when  ka  >  4  in  Fig.  4b  and  4c. 

The  middle  brand:  of  the  dispersion  curves  in  Fig.  4a  is  a 
torsional  node.  The  phase  and  group  velocities  approach  the  transverse 
wave  speed  [(1- v  )/2]  •  Due  to  axial  symmetry,  this  torsional  mode 
does  not  exist  for  the  case  of  n»0. 

The  dispersion  curves  for  the  higher  circumferential  nodes, 

n^2/3,4,.  .  .,  can  be  similarly  obtained.  Because  of  the  cubic  equation 
2 

in  ft  in  Bq.  (35) ,  there  are  three  branches  of  dispersion  curves  for 
each  value  of  n,  each  branch  corresponding  to  a  vibrational  mode. 

To  assess  the  range  of  applicability  of  Donnell's  equations  of 
notion  for  a  thin  cylindrical  shell,  the  aforementioned  dispersion 
curves  and  those  based  on  the  exact  theory  of  elasticity  for  thick- 
walled  cylinders  [11,  12]  have  been  compared.  Unlike  those  reported  in 
Ref.  11  and  12,  the  dispersion  relation  is  derived  from  three 
displacement  potentials,  and  the  calculation  of  the  dispersion  curves  is 
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programmed  for  a  hollow  cylinder  of  arbitrary  thickness.  The  results  are 
shown  as  dashed  lines  In  Figs.  3  and  4.  It  is  seen  that  a  small  error 
exists  at  the  cut-off  frequencies  for  the  higher  branches,  and  this  error 
propagates  for  ka  >  0. 

It  should  be  noted  that  based  on  the  exact  theory,  there  are  many 

more  higher  branches  which  are  not  included  in  Donnell's  theory. 

Therefore,  Donnell's  equations  should  be  applied  only  to 

analyze  tba  vibrations  of  cylindrical  thin  shells  at  a  frequency  lower 

than  the  next  highest  cut-off  frequency. 

To  include  the  effect  of  the  reacting  acoustic  pressure  for  a 
submerged  cylindrical  shell,  the  Impedance  of  8q.  (22)  should  be 
considered.  The  vanishing  of  the  total  impedance 

Z„(M)  -  Z„p(M)  ♦  Zm(k,S!)  +  ZMdc,t!)  •  0  (47) 

gives  rise  to  dispersion  curves  in  a  three  diaensional  space  formed  by  , 
real  k  and  imaginary  k.  Because  the  acoustic  impedance  of  pc  of 
water  is  30  times  smaller  than  that  of  steel,  the  effect  of  added 
impedance  *M(k,  a  )  is  expected  to  be  negligible. 

The  numerical  calculation  indicates  that  there  is  only  a  small 
difference  between  the  real  part  of  k  »  k^  +  ikj  in  Bq.  (47)  and  k 
in  Bq.  (35)  for  the  extensional  and  radial  vibrational  mode.  However, 
no  real  branches  for  the  torsional  mode  were  found  for  every 
circumferential  mode  n  >  0.  This  implies  that  the  torsional  mode  is  a 
non-radiating  mode  for  the  vibrating  submerged  cylindrical  shell 
surrounded  by  the  inviscid  fluid.  Nevertheless,  the  bending  motion  of 
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the  cylindrical  shall  caused  by  the  extansional  vibrational  mode  should 
still  contribute  to  the  radiating  energy  in  the  fluid  mediun. 


7 


NUMERICAL  RESULTS  OF  ACOUSTIC  RADIATION 


The  integrals  in  Eq.  (31)  yield  the  radiated  pressure  field  due  to 

the  harmonic  force  Fexp(-ia)  t) ,  one  for  each  circumferential  mode  n. 

The  total  radiated  pressure  is  then  obtained  by  summing  all 

circumferential  modes.  These  indefinite  integrals  are  evaluated 

numerically  for  real  wave  number  k,  and  the  summation  of  the  series  is 

truncated  when  a  desired  accuracy  is  reached. 

Of  the  three  factors  of  the  integrand  in  Eq.  (31),  1/Zn(k)  is 

regular  because  of  the  zeros  of  Z  (k)  (Eq.  (67))  are  complex  valued. 

n 

The  factor  Z  (k)  or  Z  (k)  at  r-a,  however,  has  a  singularity  at 
nr  na 

k-k„.  it  is  due  to  the  factor  (k2  -  k2)*  and  the  zero  argument  of  the 
o  o 

Hankel  functions. 

For  circumferential  mode  n  larger  than  zero,  however,  the 
singularity  at  k*ko  can  be  removed.  If  the  Hankel  functions  are 
approximated  for  small  arguments  [16]  by 

H  (z)  *  -  (i/ir)r(n)(z/2)"n  ;  H0(z)  *  (21/ir)  *n  z 

as  z  -*»  o 

Kn(z)  -v  ir(n)(z/2)‘n  ;  K0(z)  *  -  zn  z 

the  is  simplified  to 

Znr(kV'<.,V'  ioCl,o(a/n)U/r)n  ">0 
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The  integrand  of  Eq.  (31)  is  then  a  smooth  function.  For  n-0 ,  the 
singularity  is  not  removeable  because 


Zor(kVk  '  iockoa  “£lko  •  ^ 


To  minimize  the  effect  of  this  singularity  during  the  numerical 

integrations,  variable  step  Simpson's  rule  is  used. 

The  numerical  integration  in  Eq.  (31)  is  truncated  until  the  value 

of  the  indefinite  integral  converges  to  an  error  less  than  10"6.  The  infinite 

summation  over  n  is  also  truncated  when  the  sum  converges  to  error  less  than  10 

For  numerical  calculations  a  submerged  stainless  steel  cylindrical 

shell  is  selected.  The  mean  radius  of  the  shell  is  a  »  12.70  mm.  (0*5 

in.)  and  the  wall  thickness  h  ■  0.508  mm  (0.02  in.).  The 

compressional  wave  speed  of  stainless  steel  is  5.72  x  10^mm/sec; 

shear  wave  speed  c  is  3.10  *  106mn/sec,  the  Poisson  ratio  v  -  0.298 

s 

/  3 

and  its  density  Pa  ■  7.7  9/cm  •  A  receiver  is  placed  at  r  »  10a 


127.0 


(  5  in.  )  away  from  the  axis  of  the  submerged  cylindrical 


shell.  Figs.  5a,  5b,  5c  show  the  calculated  radiation  pressure  pattern 
when  the  point  force  is  driven  at  10  KHz,  50  KHz  and  100  KHz  (  ft  * 
0.139,  0.689,  1.395)  respectively. 

The  angular  distributions  of  the  pressure  fields  for  either  z-0  or 
z»r  are  plotted  in  decible  scale  relative  to  the  pressure  p  at  r  * 


10a,  9 


and  z  »  0.  The  magnitudes  of  the  pressure,  p  are 


1.25  x  10"2  dyne/ c  rr^ ,  6.14  x  10~2  dyne /cm2  and  1.31  x  io_1  dyne /on2 
at  10  KHz,  50  KHz  and  100  KHz  respectively  with  a  harmonic  driving 


force  of  F  ■  1  dyne.  For  the  two  degenerate  cases,  however,  the 


magnitude  of  the 


ring  force  is  P  ■  1/2-rra  dyne/ ca  and  that  of 


the  line  force  at  e  "0  is  P  ■  1  dyne/ca. 


\  \  \ 


.  v  \  .V'  --  *.  _*W  \-  * 


8 


DISCUSSION 


To  perform  the  numerical  integration  of  Eq.  (31)  for  the  acoustic 
radiation  over  the  real  wave  number  k  from  0  to  infinity  at  a  fixed 
frequency  and  a  fixed  circumferential  mode  n  is  equivalent  to 
integrating  along  a  horizontal  line  perpendicular  to  the  ft  axis  of  the 
dispersion  curves  in  ?igs.  3a  and  4a.  These  horizontal  lines  may 

intersect  more  than  once  on  the  real  part  of  the  complex  branches  of  the 
dispersion  curves.  The  integrand  of  2q.  (31)  evaluated  along  this 
horizontal  line  may  become  very  large  at  these  intersections,  which 

provide  a  dominant  contribution  to  the  total  pressure  intensity. 

As  shown  in  Pigs.  3a  and  4a,  the  cut-off  frequencies  of  the 

cylindrical  shell  in  vacco  increase  as  the  circumferential  mode  n 

increases.  Therefore  at  higher  frequencies  the  contribution  from 

higher  circumferential  mode  n  is  more  important  than  that  at  lower 
frequencies.  Also,  at  higher  frequencies,  the  radial  mode  is  more 
pronounced,  especially  at  small  wave  numbers.  Since  the  radiation 
pressure  depends  strongly  on  the  radial  motion  of  the  shell,  the 
pressure  may  be  higher  when  the  frequency  is  large.  The  dependence  of 
the  reference  pressure  PQ  on  the  driving  frequency  is  shown  in  Fig.  6. 

In  Pig.  5a,  the  radiation  pattern  for  a  submerged  cylindrical  shell 
driven  by  a  point  force  peaks  at  9  *  0  and  9  *  ir  »  and  dips 
approximately  at  ir  /2.  As  can  be  seen  from  the  Pourier  series  in  Sq. 
(31) ,  this  pattern  suggests  that  the  first  circumferential  mode  (n*l) 
dominates  the  vibration  of  the  shell.  In  Figure  7  the  contributions  to 
the  total  radiation  pressure  for  io  KHz  from  each  circumferential  mode 
are  shown  for  n  *  0  to  15.  The  ordinate  is  the  absolute  value  of 
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•ach  term  in  tha  series  of  Eq.  (31)  for  r  ■  10a(127  ram) ,  9  *  0°,  z  - 
0  or  10a.  The  n  ■  1  mode  (translational  mode)  is  indeed  dominant. 

Fig.  5b  exhibits  two  peak  pressures  at  9  ■  0  and  0  ■  2  ir/3  for 
a  submerged  cylindrical  shell  driven  by  a  point  force  at  50  KHz. 
Because  of  the  cos  n€  dependence  of  the  radiated  field,  these  two  peaks 
seem  to  be  the  effect  of  a  term  with  n  ■  3/2.  In  terms  of  the  integral 
circumferential  mode  number,  the  combination  of  n  ■  1  and  n  »  2  mode 
could  be  the  dominant  source  for  the  radiation  field.  This  is 
substantiated  by  the  contributions  from  each  mode  as  shown  in  Fig.  7b. 

At  100  KHz  in  Fig.  5c,  there  are  no  distinct  peaks  and  dips. 
Therefore  not  a  single  or  two  circumferential  modes  are  dominant.  This 
is  verified  by  the  results  shown  in  Fig.  7c. 

As  for  the  submerged  cylindrical  shell  driven  by  a  ring  source  at  z 
■  0,  because  of  the  axi symmetrical  source,  the  radiation  pattern  is 
independent  of  6  as  shown  in  the  lower  half  of  Figs.  5a,  5b,  and  5c. 
Similarly  for  the  shell  driven  by  a  line  source  at  9  *  0,  the  radiation 
field  is  independent  of  the  z-coordinate. 

To  summarize,  the  radiation  patterns  of  the  pressure  field  in  a 
fluid  surrounding  a  circular  cylindrical  shell  are  strongly  influenced 
by  the  nature  of  the  force  that  excites  the  shell  (  a  point  force,  a 
ring  force,  and  a  line  force),  and  are  dependent  on  the  frequencies  of 
the  driving  force.  At  lower  frequencies,  the  circumferential  modes  n  * 
1  and  2  are  dominant  and  the  angular  distribution  of  the  pressure  peaks 
at  9  ■  0  and  *  and  dips  near  9  *  tt/2  for  n  ■  1.  At  higher 
frequencies,  the  lowest  four  or  five  circumferential  modes  are  equally 
important  in  the  contribution  to  the  total  wave  field.  Since  the 


highest  branch  of  the  dispersion  curve  of  the  shell  based  on  Donnell's 
equation  is  predominantly  radial  motion  (at  lower  wave  numbers) ,  and  the 
radiation  field  is  excited  mainly  by  the  radial  motion  of  the  shell ,  the 
pressure  field  is  found  to  be  stronger  for  a  driving  force  with  higher 
frequencies. 
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Fig.  2.  Relation  between  the  wave  number  kQ  and  k. 

Fig.  3.  Dispersion  curves  of  extensional  waves  in  a  thin  shell  (n  *  0) 

(a)  Frequency  versus  wave  number 

(b)  Phase  velocity  versus  wave  number 

(c)  Group  velocity  versus  wave  number 


Fig.  4.  Dispersion  of  flexural  waves  in  a  thin  shell  (n  *  1) 

(a)  Frequency  versus  wave  number 

(b)  Phase  velocity  versus  wave  number 

(c)  Group  velocity  versus  wave  number 

Fig.  5.  Angular  distribution  of  radiation  pressure  at  r  *  10  a  for  a  submerged 
steel  shell  of  mean  radius  a  *  12.7  mm.  The  shell  is  driven  by  a 

point  force  of  one  dyne  at  r  *  a,  a=*  0,  z  *  0  (the  upper  half  circle). 

Solid  lines  are  for  pressure  at  z  *  0  and  dashed  lines  for  pressure  at 
z  »  10  a.  Pressure  fields  generated  by  a  line  force  at  r  *  a,  9  *  0 
are  shown  in  the  lower  half  circle  by  dotted  lines.  All  pressures  are 
measured  in  db  relative  to  the  peak  reference  pressure  pQ  at  r  =  10  a, 

Q  *  0,  z  *  0. 

(a)  Shell  driven  at  10  KHz,  pQ  =  0.00125  dyne/cm2. 

(b)  Shell  driven  at  50  KHz,  pQ  =  0.00614  dyne/cm2. 

(c)  Shell  driven  at  100  KHz,  pQ  =  0.0131  dyne/cm2. 


Variation  of  the  reference  pressure  with  the  driving  frequency. 

Absolute  values  of  the  pressure  field  at  (10a,  0,  0)  in  solid  line 
or  at  (10a,  0,  10a)  in  dashed  lines  due  to  each  component  of  the 
circumferential  mode  (n  *  0  to  15). 

(a)  Shell  driven  at  10  KHz. 

(b)  Shell  driven  at  50  KHz. 

(c)  Shell  driven  at  100  KHz. 
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Fig.  4(c)  Dispersion  of  flexural  waves  in  a  thin  shell 
Group  velocity  versus  wave  number 


Angular  distribution  of  radiation  pressure  at  r  3  10  a  for  a  submerged 

steel  shell  of  mean  radius  a  3  12.7  mm.  The  shell  is  driven  by  a 

point  force  of  one  dyne  at  r  =  a,  3  »  0,  z  =  0  (the  upper  half  circle) 

Solid  lines  are  for  pressure  at  z  =  0  and  dashed  lines  for  pressure  at 

z  3  10  a.  Pressure  fields  generated  by  a  line  force  at  r  =  a,  e  3  0 

are  shown  in  the  lower  half  circle  by  dotted  lines.  All  pressures  are 

measured  in  db  relative  to  the  peak  reference  oressur*  o  at  r  =  10  a 

o  ’ 

3  s  0,  z  3  0. 

(a).  Shell,  driven  at  10  KHz,  p_  3  0.00125  dyne/cn2. 


Fig.  5.  Angular  distribution  of  radiation  pressure  at  r  *  10  a  for  a  submerged 

steel  shell  of  mean  radius  a  *  12.7  mm.  The  shell  is  driven  by  a 

point  force  of  one  dyne  at  r  =  a,  9  =  0,  z  *  0  (the  upper  half  circle) 

Solid  lines  are  for  pressure  at  z  *  0  and  dashed  lines  for  pressure  at 

z  *  10  a.  Pressure  fields  generated  by  a  line  force  at  r  =  a,  3  =  0 

are  shown  in  the  lower  half  circle  by  dotted  lines.  All  pressures  are 

measured  in  db  relative  to  the  ceak  reference  pressure  p  a*  r  =  10  a 

o  ’ 

9  *  0,  z  *  0. 

(b)  Shell  driven  at  50  KHz,  pQ  =  0.00614  dyne/cm^ . 
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Angular  distribution  of  radiation  pressure  a.t  r  =  10  a  for  a  submerged 
steel  shell  of  mean  radius  a  =  12.7  mm.  The  shell  is  driven  by  a 
point  force  of  one  dyne  at  r  =  a,  9  a  0,  z  -  0  (the  upper  half  circle) 
Solid  lines  are  for  pressure  at  z  =  0  and  dashed  lines  for  pressure  at 
z  a  10  a.  Pressure  fields  generated  by  a  line  force  at  r  =  a,  a  *  0 
are  shown  in  the  lower  half  circle  by  dotted  lines.  All  pressures  are 
measured  in  db  relative  to  the  peak  reference  oressure  pQ  at  r  =  10  a, 
a  a  0,  z  a  0. 

(c)  Shell . driven i.atjOO  KHz_.  o  *  0.0131  dvne/r-r 


Absolute  values  of  the  pressure  field  at  (10a,  0,  0)  in  solid  line 
or  at  (10a,  0,  10a)  in  dashed  lines  due  to  each  component  of  the 
circumferential  mode  (n  =  0  to  15). 

(a)  Shell  driven  at  10  KHz. 


n 


Absolute  values  of  the  pressure  field  at  (10a,  0,  0)  in  solid  line, 
or  at  (10a,  0,  1 0a )  in  dashed  lines  due  to  each  component  of  the 
circumferential  mode  (n  *  0  to  15). 


Absolute  values  of  the  pressure  field  at  (10a,  0,  0)  in  solid  line 
or  at  (10a,  0,  10a)  in  dashed  lines  due  to  each  component  of  the 
circumferential  mode  (n  *  0  to  15). 

(c)  Shell  driven  at  100  KHz. 
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Abstract 

A  submerged  shell  of  finite  length  is  driven  by  a  concentrated 
force  in  steady  state,  and  radiates  acoustic  waves  into  the  surrounding 
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expressed  in  terms  of  Fourier  series,  and  the  unknown  coefficients  are 
shown  to  satisfy  a  system  of  integral  equations.  These  equations  can  be 
reduced  to  a  system  of  algebraic  equations  and  solved  numerically. 
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1.  INTRODUCTION 

The  acoustic  radiation  of  a  submerged  vibrator  is  a 
subject  of  many  studied  [4,9,12,15],  The  importance  of  this 
topic  lies  in  its  application  to  the  noise  reduction  of  a 
vibrator  and  to  the  understanding  of  the  dynamic  response  of 
submerged  vehicles.  The  common  approach  of  analyses  can  be 
divided  into  two  categories;  one  is  the  application  of  Green's 
function  integrated  over  the  entire  surface  of  the  submerged 
object  for aF  prescri bed  excitation  at  the  fluid-vibrator  inter¬ 
face  [4],  the  other  is  to  solve  the  coupled  equations  for  the 
vibrator  and  the  surrounding  fluid  medium  [12],  The  first 
approach  eliminates  the  complexity  of  the  interaction  between 
the  reactive  pressure  and  the  vibrating  object.  The  second 
is  usually  very  complicated  except  for  those  vibrators  with 
simple  geometries  such  as  a  cylindrical  or  spherical  shell. 

In  this  paper,  a  theoretical  study  of  the  acoustic  radia¬ 
tion  from  a  submerged  finite  cylindrical  thin  shell  driven  by 
a  point  harmonic  force  is  presented.  Donnell's  equations  of 
motion  for  the  circular  cylindrical  thin  shell  are  adopted  to 
analyze  the  deformation  of  the  thin  shell,  and  the  scalar  wave 
equation  for  the  surrounding  fluid  medium  is  used  to  investigate 
the  acoustic  radiation  field  induced  by  the  vibration  of  the  shell. 


The  cylindrical  shell  of  finite  length  is  sealed  by  two 
flat  end  caps.  For  simplicity,  the  deformation  of  the  flat 
end  caps  are  neglected  in  the  analysis.  Thus  the  end  caps 
move  like  a  rigid  body  synchronous  to  the  longitudinal  motion 
of  the  thin  shell.  The  piston-like  motion  generates  additional 
acoustic  waves  into  the  surrounding  medium.  The  induced  acoustic 
pressure  field  at  the  fluid-shell  interface  resulted  in  a  react¬ 
ing  force  for  the  vibration  of  the  cylindrical  thin  shell. 

This  interaction  between  the  motion  of  the  thin  shell  and  the 
acoustic  pressure  renders  the  analysis  very  complex. 

In  the  first  part  of  this  paper,  Donnell's  equations  of 
motion  for  the  cylindrical  thin  shell  are  introduced.  Because 
of  the  symmetric  condition  of  the  point  excitation  in  the 
circumferential  direction,  the  displacement  of  the  shell  can 
be  expressed  in  the  form  of  Fourier  series  in  the  angular 
coordinate.  The  pressure  wave  field  is  expanded  in  terms  of 
the  spherical  harmonic  as  it  forms  a  spherical  wave  front  in 
the  far  field.  In  the  next  section,  the  scalar  wave  equation 
for  the  pressure  field  is  solved  by  using  the  method  of  Green's 
function  for  the  Helmholtz  equation.  Through  the  continuity 
condition  of  the  fluid  particle  velocity  and  the  normal  veloci¬ 
ties  of  the  thin  shell,  the  pressure  field  is  expressed  in 
terms  of  the  deformation  of  the  shell.  In  the  third  section, 
the  appropriated  boundary  conditions  for  the  finite  cylindrical 
thin  shell  are  applied.  In  the  fourth  section,  the  coupled 
equations  between  the  motion  of  the  shell  and  the  fluid  medium 


are  reduced  to  an  inhomogeneous  integral  equation  for  the  dis¬ 
placements  with  a  separable  kernel.  This  integral  equation 
is  solved  by  reduction  to  a  system  of  algebraic  equations. 

Then  the  radiated  pressure  is  readily  obtained  from  the  displace 
ment  of  the  cylindrical  thin  shell.  In  the  last  section,  the 
difference  between  the  continuous  and  discrete  wave  number  is 
discussed.  In  addition,  some  possible  modifications  for  the 
analysis  of  a  submerged  finite  cylindrical  thin  shell  with 
different  end  caps  are  also  discussed. 
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2.  MOTION  OF  A  SUBMERGED  FINITE  CYLINDRICAL  THIN  SHELL 

Consider  a  cylindrical  thin  shell  of  finite  length  2L 
sealed  with  two  flat,  rigid  caps  immersed  in  an  infinitely 
extended  fluid  medium.  The  cylindrical  thin  shell  is  excited 
by  a  harmonic  point  force  F  exp(-icot),  where  F  is  the  ampli¬ 
tude  of  the  driving  force,  and  to  is  the  circular  frequency  of 
the  force  acting  in  the  radial  direction  against  the  inner 
surface  of  the  shell.  The  two  sealed  end  caps  are  assumed 
to  be  weightless  thin  plate  with  negligible  rigidity.  For 
simplicity,  the  deformation  of  the  end  caps  are  neglected  in 
the  following  analysis.  Thus  the  boundary  conditions  for  the 
vibrating  cylindrical  thin  shell  are  assumed  free  at  both  ends. 
It  is  assumed  that  the  cylindrical  shell  is  in  steady  state 
vibration  and  a  time  dependence  of  the  form  exp(-iwt)  for  all 
dependent  variables  is  assumed  throughout  the  derivation. 

Because  of  the  circular  cylindrical  geometry,  Donnell's 
equation  for  cylindrical  thin  shell  is  applied  to  analyze  the 
dynamic  behavior  of  the  shell.  This  equation  of  motion  had 
been  studied  extensively  [6,7,9,14].  The  omission  of  higher 
overtones  of  the  radial,  torsional  and  longitudinal  vibrational 
modes  of  cylindrical  shell  in  Donnell's  theory  limits  its 
applicability  below  the  cut-off  frequency  of  its  lowest  simple 
thickness-shear  vibrational  mode  [3],  Thus,  the  present 
derivation  is  valid  only  at  lower  frequencies. 

The  vibration  of  the  cylindrical  shell  radiates  acoustic 


waves  into  the  surrounding  fluid  medium.  Because  of  the 


finite  dimension  of  the  vibrator,  the  radiation  is  in  the 
form  of  spherical  waves  in  the  far  field.  The  radiated  pressure 
wave  is  thus  expressed  in  the  spherical  coordinate  system.  The 
addition  theorem  is  used  to  transform  the  spherical  wave  func¬ 
tions  into  those  in  the  cylindrical  coordinates  in  order  to 
match  the  wave  field  with  the  vibration  of  the  cylindrical 
shell  . 

For  a  thin  cylindrical  shell,  as  shown  in  Fig.  1,  the 
equations  of  motion  are  [8] 
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where  r,0,z  are  the  cylindrical  coordinates,  u,v,w  are  the 
displacements  in  the  longitudinal,  tangential  and  radial 
directions,  respective,  v  is  the  Poisson  ratio,  E  is  the 


Young's  modulus,  Cp  =  [E/p $ ( 1  -v2 ) is  the  compressi onal  wave 


speed  of  the  plate,  pg  is  the  mass  density,  h  is  the  thickness. 


and  a  is  the  mean  radius  of  the  shell. 


On  the  right-hand  side  of  Eq.  (la),  the  acoustic  pressure 
fj(a,0,+L,t)  and  p(a,e,-L,t)  at  two  end  caps  form  two  reacting 
longitudinal  forces.  On  the  right-hand  side  of  the  third 
equation,  the  first  term  is  due  to  the  radially  driven  point 
force  and  the  second  is  the  reacting  acoustic  pressure  on  the 
cylindrical  fluid-shell  interface. 

Due  to  the  symmetry  of  the  driving  point  force  F5(e)/a, 
the  displacements  u,  v,  w  and  the  Kronecker  delta  function 
6(0)  can  be  expressed  as 
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and 
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where  the  index  m  for  the  circumferential  vibrational  mode 
indicates  the  number  in  pairs  of  the  nodal  points  along  the 
circumference  of  the  circular  cylindrical  shell  (see  Fig.  2). 
Because  of  the  spherical  wave  front  of  the  radiation  at 


far  field,  the  pressure  wave  P(r,e,z)  is  conveniently  expressed 
in  terms  of  spherical  harmonics 


p(r,e,z)  =  p(p,e,<|>)  =  l  l  A  Pnm(cos  4> )  hn  '  ( kp  )cos  me 

n=0  m=0 


where  p,e,<j>  are  the  spherical  coordinates,  Pnm(cos  <t>)  is  the 
Legendre  function,  h  ^(kp)  is  the  spherical  Hankel  function 
of  the  first  kind  [1],  Anm  are  the  unknown  coefficients  and 
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Substituting  Eqs.  (2),  (3)  and  (4)  into  Eq.  (1)  and  utilizing 
the  orthogonality  properties  of  the  sinusoidal  functions,  we 
have 
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where  H(z+L),  H(z-L)  in  Eq.  (6c)  are  the  Heaviside  step  func¬ 
tions  and  the  equivalent  reacting  pressures  in  Eq.  (6a)  are 
obtained  by  integration  over  a  surface  area  substended  by  a 
small  angle  A6  on  both  end  caps. 

To  solve  the  simultaneous  equations  of  motion  (6),  the 
axial  dependence  of  the  displacements  in  Eqs.  (2)  are  assumed 
to  be 


*  /:  y*)eUzdx 

V*)  *  /:  Vm(x)eU2dx  (7) 

W  ( z )  =  r  Wm(x)eiXzdx 
m  m 

where  x  is  the  wave  number  in  the  axial  direction.  The  coef¬ 
ficients  Um(x),  VM(x)  and  Wm ( X )  in  Eq.  (7)  can  be  obtained  by 
using  the  generalized  Fourier  transform 

y*>  ■  k/:!x.<2>e',xz” 

».<*>  *  <8) 


Substituting  Eq.  (7)  into  Eq.  (6)  and  utilizing  the  generalized 
Fourier  transform  (8),  we  have  the  Donnell's  equations  (1) 
expressed  in  the  wave  number  domain  (m,x) 
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Replacing  the  expression  [  - . . . J in  the  r-{9ht  hand  side 
Eq.  (9)  by  the  Integral  defined  in  Eq.  (7) 
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3.  ACOUSTIC  WAVES  IN  THE  FLUID  MEDIUM 

The  radial  motion  of  the  vibrating  cylindrical  shell 
associated  with  the  longitudinal  motion  of  the  end  caps  gen¬ 
erate  a  radiating  pressure  field  in  the  surrounding  fluid 
medium.  If  this  disturbance  in  the  inviscid  fluid  medium  is 
infestimal,  the  pressure  field  P(p,e,*,t)  satisfies  the 
linear  wave  equation  [2] 

v2p(P,e,<j,,t)  =  J-3-p(P ,e,*,t)  ,  cu  *  /nr  (12) 

w  yj  WWW 

where  p,0,<j>  are  the  spherical  coordinates,  v2  is  the  Laplacian 
operator  and  Cw»  pw,  xw  are  the  wave  speed,  mass  density  and 
Lame's  constant  of  the  fluid  medium,  respectively.  In  steady 
state 


P(p  >e  »<fr  ,t)  =  p(P  ,0,*)e~iwt  ,  (13) 

Eq.  (12)  is  transformed  into  Helmhortz  equation 

v2p(p,0,$)  +  k2p(p  ,9 , 4» )  =0  ,  k  =  «o/Cw  (14) 

subjected  to  the  boundary  conditions  at  the  fluid-shell  inter¬ 
faces  where  continuity  requires  the  shell  and  fluid  particle 
velocities  to  be  equal. 


Equation  (14)  is  solved  by  the  method  of  Green's  function 
so  as  to  include  the  yet  undetermined  boundary  conditions  at  the 


shell  surface.  Green's  function  satisfies  an  inhomogeneous 
three  dimensional  equation. 


?2G(r|r0)  +  k2G(r|r0)  =  <$(r  -  rQ) 


In  an  unbounded  medium  [10],  the  solution  is 
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where 


R2  *  r2  +  rQ2  -  2rrQ  cos(e-0o)  +  (z  -  zQ)2 

=  p2sin2i>  +  P02sin2$0  -  2ppQsin  sin  <j>0  cos(e-0Q) 


+  (p  cos  <j>-  pQ  cos  <f>Q)2 


(17) 


is  the  distance  between  the  source  point  (pq  »90  [  ( r0  »90  *z0) 

in  the  cylindrical  coordinate]  and  the  observation  point 
(p .0 »+)[(r ,0 ,z)  in  the  cylindrical  coordinate].  The  Green's 
function  G(rjrQ)  of  Eq.  (16)  can  be  expanded  in  terms  of  the 
spherical  harmonics  [11] 


i  JkR 
G(r|r0)  -  ji  K~ 


ik.  (1) 
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( kR) 


r1  l  ( 2n+i )  i  em|^a 
4"n=0  m«0 


( n-m 


P„m(cos  "‘(cos  ♦)( 


m. 


!Jn(kro)hn<1)(kr)>  r  -  r( 


(:.(kr)h„<1)(kr„),  r 


<  r. 


(18) 


To  obtain  the  solution  of  Eq.  (14),  we  multiply  Eq.  (14) 

by  6(r| r  )  and  Eq.  (15)  by  p(r)  and  subtract,  exchanging  r 

for  r  at  the  same  time 
-o 


G(r  | 


r  )  v  2  p  ( r  )  - 

-o'  o  K ' -o' 


p<r0>v0- 


G  ( r 


r0)  - 


(r0)s(r0-r)  •  09) 


Integrating  this  equation  over  all  the  source  coordinate  r  , 
we  obtain 


(p(r)  r  outside 

///  [S(r|r  )v  2p(r  )-p(r  )v  2G(r|r  )]dV  =  (  ~  the  shell 

v  (o  if  r  inside 

0  (  theshel 1 

(20) 

Then  tne  Green's  Theorem  is  used  to  simplify  the  volume  integral 
of  Eq.  (20)  into  a  surface  integral 


#s  C-fitrlEoJV^o) 


+  P(r0)v0S(r|r0)]-dA 


(p(r)  r  outside 
=  (  ~  the  shell 

(o  r  inside  the 
(  shell 


(21) 

The  gradient  of  the  surface  element  dAQis  pointing  outward  from 
the  cylindrical  shell  into  the  infinitely  extended  fluid  medium. 
Since  the  fluid-shell  interface  are  made  of  three  distinct 


surfaces  of  the  vibrating  body,  one  cylindrical  shell  and  two 
flat  end  caps,  the  total  surface  area  sQ  of  Eq.  (21)  is  then 
subdivided  as 


»<:)  ■  U. l/0  Vd9odzo  +  Mo  -zrodeodro  +  /0/o  (-?2)rodeodro:i 
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0 
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0 

,  r^ikR  a p ( r ^ )  .  r«ikRy-r 

+  Vo  0K - 5^7 - p(ro >9l^  (t-)12  =_Lrod9odro  (22) 


The  pressure  p  is  related  to  the  fluid  particle  velocities 
vf  by  the  balance  of  momentum  [2] 


-vp  *  pw-f 


(23) 


and  the  fluid  particle  velocities  are  equal  to  the  normal  com¬ 
ponents  of  the  velocities  of  the  shell  at  the  fluid-shell 
interface  s^ 


t*f]r  .a,n=e  -  *<a-vV 

0  -  -  r 


Cvf]z-L,n-ez  =  °*a,Bo,zo' 


(24) 
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The  gradient  of  pressure  VQp(r0)  in  Eq.  (22)  is  then  written  as 
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Thus,  Eq.  (22)  is 
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(25) 


(26) 


Substituting  the  expressions  in  Eqs.  (2),  (4),  (7)  and  (18) 
i nto  Eq .  (26) ,  we  have 
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Note  that  Eq.  (27)  can  be  rewritten  as 
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by  using  the  orthogonality  properties  of  the  spherical 
Then  in  the  matrix  form,  Eq .  (30)  is 


[D]{A>  =  /“oo[C(x)]{U}dX  for  fixed  m 

where 
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and  {A},  (U)  are  defined  in  Eq.  (11). 

Due  to  the  feet  that  Z™nI  =  z”„t  =  z”nl  =  Z4nm  = 

=  Z,„m  =  0  when  m  >  n  or  m  >  £ ,  the  matrices  A,  C  and 


harmoni cs 

(31) 
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(32b) 
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4.  BOUNDARY  CONDITIONS 

The  boundary  conditions  for  the  present  finite  cylindrical 
thin  shell  are  assumed  to  be  two  free  ends  provided  this  shell 
is  sealed  by  two  end  caps  with  negligible  mass.  For  cylindrical 
shell,  the  free  ends  implies  [13] 


longitudinal  force:  N. 


Eh  fau 


-V2  3Z 


3  V  Wf[ 

a 30  ~  aJJ 
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shear  force:  Sz  -  |f] 


0(1  — 


at  z  =  ±L 


^  ♦  0e]=° 
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bending  moment:  M,  *  -D  ft  +  tt 
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torque:  T  *  -D~-y  *  TT 
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where  D  *  Eh3/[1 2( 1 -v2 ) ]  is  the  flexural  rigidity  of  the  cylind¬ 
rical  shell.  After  replacing  the  di spl acements  by  their  Fourier 
expansion  defined  in  Eq.  (2)  and  utilizing  the  orthogonality 
property  of  the  sinusoidal  function,  the  above  boundary  condi¬ 
tions  can  be  simplified  as 
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(35b) 
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Substituting  the  above  boundary  conditions  (35)  into  the 
corresponding  expression  on  right-hand  side  of 

Eq.  (9);  e.g.  dUm/dz  in  Eq.  (35a)  into  Eq .  (9a),  d"Vm/ dz  in 
Eq.  (35b)  into  Eq.  (9b),  d2Wm/dz2  and  d3Wm/dz3  in  Eq.  (35c) 
and  (35d),  respectively,  into  Eq.  (9c);  the  integral  equation 
of  Eq.  (10c)  should  be  rewritten  as 


[S]{U>  =  [Q]  {A}  +  {B}  +  ^r/IC0[K’(X,5)]{U(5)}sin[(x-c)L]d5  (36) 

where  the  kernel  [K'(x,0]  is 
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5.  DISPLACEMENTS  OF  THE  CYLINDRICAL  THIN  SHELL 

Substituting  Eq.  (31)  into  Eq.  (36),  we  have  the  integral 
equation  for  the  displacements  {U} 

[S(X)]<U(X)>  ■  CQ(X)] CD]_1/"„ CC(x)]{U(c)}d«  +  {B ( X ) } 

+  |^/ra>[K'(x,5)]{U(5)}sin[(x-c)L]dc  (38a) 

or 

[S(x)]{U(x)>  =  (B(x ) }  +  /_„[ K(x,0]{UU)}<U  (38b) 

where 

[K(X.C)]  -  J7Sil>[(x-5)L][K'(X.E)]  ♦  [Q(X)][D]‘1[C(C)]  .  (38c) 

To  solve  for  the  displacements  {U},  we  notice  that  this  integral 
equation  (38b)  has  a  degenerate  kernel 

[K(x.O]  *  ^s1nC(x-OL]([q'(x)]  +  [ICJU)]  +  xHKJU)] 

+  [Q(x)][0]'1[C(5)]  (39a) 

where 
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[KJ(X)] 


(39b) 


|o 


Thus,  this  integral  equation  can  be  solved  by  reduction  to  a 
system  of  algebraic  equations  [10]. 

Substituting  (39a)  into  Eq.  (38b),  we  have  the  following 
system  of  integral  equation  for  this  displacement  { U } 
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TU(  A) }  -  l  [6i(x)]/c#  [Ki(5)]{U(5)}d5=  (B'(x)J  ( 

1*1  -  00 

where 

[G,(x)]  ■  ^r[S(x)]-1[K1»(x)]s1n(xL)  ( 
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[Gj(x)3  *  S ( x ) ]  ^sln(xL)  ( 

[G.(x)l  ■  ^■rs(x)]-1cos(xL)  ( 


[6 5(x)]  =  f7[S(x)]-,sin(xL) 
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[K-|  (5)]  =  [ I ]cos ( ?L ) 

[K2U)]  -  [I]sin(§L) 

C  K3 ( 5 ) ]  =  [K”(5 ) ]cos (?L) 

[K4U)]  =  [K"<^)]sin(5L) 

U5U)]  *  [Kjj($)]cos(gL) 

C Kg ( e ) 3  *  [K|(g)]s1n(gL) 

[k7(5)]  -  CC(e)] 

{ B 1 ( X ) }  =  [S(X)]-1{B(X)} 

If  we  define 

<ei>  s  Jl«CKiU)]{UU)>di;  ,  1-1  ,2,. ..7 


the  Integral  equation  (40)  Is  then  rewritten  as 


(43) 


{U(X)}  =  X  [  G  .  ( X )  ]  { e  . }  +  { B  *  (  X )  > 
i  =1  1  1 

The  column  matrix  {e^}  in  Eq.  (42)  are  unknown  constants, 
since  the  function  {U(x)}  is  unknown.  To  solve  for  the  unknown 
constants  {e..},  we  substitute  Eq.  (43)  into  the  integral  equa¬ 
tion  (40)  and  obtain  a  linear  combination  of  the  functions 
[  6 1  ( X )  ] 

l  CS1(x)]{{ei>  -/.JXU)]  C  B  •  (e)  >  +  \  [Gk(c)]{ek}]d5}  =  0  . 

i -1  1  11  k  =  l  K  K 

(44) 

Because  of  the  linear  independence  of  the  functions  [G.(x)] 
within  the  interval  {X;  -•»  <  X  <  «},  the  above  Eq.  (44)  implies 

7 

{ei>  -/i.CK1(5)][(B,(5)}  +  l  [GkU)]{ek}]dc  =  0  , 

k  —  1 

i  -  1  .2 . 7.  (45) 

This  introduces  two  calculable  integrals  which  may  be  denoted 
as  follows 

C^j]  -7!.CK1(5)]C6J(e)]dc  (46) 

[Bi  3  =  [Ki (s):TB ' (e) }de  .  (47) 

Then,  Eq.  (45)  can  be  rewritten  as 


0 


(48) 


{ei}  "  k^Caik]{ek}  '  {6i}  = 

or 
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l  Cl  -  «ik^ek>  =  { B-j }  »  1  3  1  ,2,. ...7  (49) 

k-1 

where  [I]  is  a  unit  matrix  of  order  3  and  [aikl  is  also  a  matrix 
defined  in  Eq.  (46).  Note  that  we  have  obtained  a  system  of 
linear  algebraic  equations  for  the  determination  of  the  unknown 
constants  {ek>.  Solving  these  simultaneous  Eqs.  (49)  with  a 
21  x  21  coefficient  matrix  [I  -  a.k],  the  solution  of  the  dis¬ 
placements  { U }  are  then  given  by  Eq.  (43). 

We  note  that  the  determinant  of  the  coefficient  matrix 


E(oj)  =  [I  -  aikU)] 


I  1  1  (tu) 

**a-|  2  ( u ) 

-t*l  7(ai) 

— ex 2 i  (“) 

1 -a 22^“ )  *  *  * 

-a27(ai) 

(50) 


U) 


is  a  transcendental  function  of  the  frequency  <d.  Therefore,  it 
Is  possible  to  determine  numerically  the  resonant  frequencies 
of  the  submerged  finite  cylindrical  thin  shell  such  that 
det[E(w)]  *  0.  In  such  case,  the  system  of  Eq.  (49)  does  not 
have  a  unique  solution.  In  another  words,  at  least  one  of  the 
{ek)'s  should  be  assigned  arbitrary,  and  the  remaining  {e.,}'s 
can  be  determined  accordingly. 


For  those  frequencies  de t[ E ( oj ) ]  f  0,  we  may  determine 
uniquely  the  unknown  constants  {e^}  according  to  the  frequency 
and  the  position  of  the  excitation  specified  in  Eq.  (10). 
Substituting  these  constants  { e k }  into  Eq.  (43),  we  obtain 
the  displacements  of  the  submerged  finite  cylindrical  thin 
shell  {U(x)}  in  the  wave  number  domain  x.  Then,  the  inverse 
Fourier  transform  (7)  and  the  summation  of  Fourier  series  (2) 
are  performed  to  obtain  the  displacements  {u(a,e,z)}  in  the 
spatial  domain. 


6.  ACOUSTIC  RADIATION  FIELD  OF  A  SUBMERGED  FINITE  CYLINDRICAL 


SHELL 

Once  the  displacements  {U(x)>  in  the  wave  number  domain  x 
are  determined  by  Eq.  (43),  the  expansion  coefficients  Anm  for 
the  acoustic  radiation  field  can  be  obtained  from  Eq.  (31) 

{A}  =  [D]’1/00  [C(x)]{U(x)}dx  for  fixed  m  .  (51) 

*  «•  oo 


Substituting  these  coefficients  Anm  into  Eq.  (14),  the  acoustic 
radiation  field  outside  the  vibrating  cylindrical  shell  is 
obtai ned 


Pfr.e.z.t)  -  P(p,..*.t)  =  \  j  AnmPj  '(7frf!r]h‘1)(k/r*«*) 

n=o  m=u 


cos  mee 


•iwt 


(52) 


7.  DISCUSSION 


It  is  known  that  the  vibration  of  an  infinitely  long 
cylindrical  shell  may  be  described  by  a  continuous  wave  number 
in  the  axial  direction.  These  broadband  spatial  spectra  degen¬ 
erate  to  a  discrete  set  of  wave  numbers  when  the  dimension  of 
the  vibrator  is  finite.  In  the  previous  derivation,  the 
displacement  { U ( x ) >  in  Eq.  (7)  and  the  acoustic  field  P(r) 
in  Eqs.  (51)  and  (52)  indicate  a  continuous  wave  number  x  should 
be  included.  These  facts  are  mainly  attributed  to  the  reacting 
pressure  field  on  the  two  flat  end  caps.  These  vibrating  end 
caps- generate  the  acoustic  pressure  wave  propagating  in  the 
infinitely  extended  fluid  medium.  These  acoustic  waves  will 
form  a  spherical  wave  front  in  the  far  field  and  have  wave 
number  with  Infinite  kinds  of  projection  in  the  axial  direction. 
Thus  the  wave  number  x  should  indeed  have  a  continuous  set. 
However,  the  previous  study  on  the  acoustic  radiation  from  an 
Infinitely  long  submerged  cylindrical  shell  [9]  indicates  the 
fluid  medium  acts  like  a  resistive  impedance  when  the  waves 
in  the  cylindrical  shell  have  the  wave  number  in  the  axial 
direction  larger  than  the  wave  number  in  the  fluid  medium 
u/c  (x  >  w/c  ) .  Thus,  the  dominant  contribution  of  the  in- 
finitely  extended  integration  of  X  for  the  displacements  { U } 
in  Eq.  (7)  and  the  radiated  acoustic  field  in  Eq.  (51)  should 
emerge  from  the  integration  when  |x|  <  oj/c  .  Then,  the  infinitely 

W 

extended  integration  can  be  truncated  in  according  to  the  desired 
accuracy.  Similarly,  the  infinite  series  of  the  circumferential 


mode  in  Eqs.  (2)  and  (52)  may  also  be  truncated  because  the 
increasing  stiffness  of  the  cylindrical  thin  shell  at  the 
higher  circumferential  vibrational  modes. 

In  the  present  analysis,  two  flat  end  caps  for  the  finite 
cylindrical  shell  have  been  assumed.  By  neglecting  the  defor¬ 
mation  of  the  flat  end  caps,  it  is  reasonable  to  assume  the 
boundary  conditions  for  this  finite  cylindrical  shell  to  be 
two  free  ends.  Should  the  end  caps  be  replaced  by  two  elastic 
flat  plates  or  two  spherical  caps,  the  boundary  conditions 
should  be  modified  to  the  continuity  of  the  displacements  and 
the  stresses  at  the  junctions  between  the  finite  cylindrical 
shell  and  the  sealed  caps. 
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Abstract 

The  phase  velocities  and  attenuation  constants  of  eighteen 
polymeric  materials  are  measured  by  ultrasonic  methods.  Both 
longitudinal  (pressure)  waves  and  transverse  (shear)  waves  are 
Investigated.  The  attenuation  constants  are  found  to  be  approximated  by 
a  quadratic  polynomial  of  the  frequency  variable  within  the  range  of  1 
to  10  MHz. 
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ATTENUATION  CONSTANTS  AND  PHASE  VELOCITIES 
OF  EIGHTEEN  POLYMERIC  MATERIALS 


C.P.  Chen 


Experiments  were  conducted  to  measure  the  frequency  dependent  attenuation 
constants  and  phase  velocities  for  both  longitudinal  and  shear  elastic  waves 
in  eighteen  polymeric  materials.  The  broadband  pulse-echo  method  was  used  for 
all  measurements.  Frequency  range  was  from  1  to  10  MHz. 

Ultrasonic  piezoelectric  transducers  are  used  as  transmitter  and  receiver. 
The  schematic  diagram  of  the  instruments  is  shown  in  Figure  1,  and  a  typical 
record  of  the  received  signals  is  shown  in  Figure  2.  Attenuation  measurements, 
the  first  two  echos  of  ultrasonic  signal  propagating  back  and  forth  in  the 
specimen  are  digitized  and  windowed  by  the  digital  sampler.  The  sampled  data 
are  then  analyzed  by  applying  a  FFT  routine  to  obtain  the  Fourier  amplitude 
spectra  of  two  echos.  The  attenuation  is  calculated  by  the  following  formula: 


,fV  NP  _  1  , 

a(f)  cm  '  2d  logi 


XJJT 


where  A-j(f)  is  the  Fourier  amplitude  spectrum  of  the  first  echo  and  A2(f) 
is  the  Fourier  amplitude  spectrum  of  the  second  echo,  both  at  the  same  fre¬ 
quency  f,  and  d  is  the  thickness  of  the  specimen. 

For  phase  velocity  measurement,  phase  spectrum  analysis  technique  is  used 
to  calculate  the  frequency  dependence  of  phase  velocity.  The  procedures  of 
measurement  is  described  in  the  paper  by  W.  Sachse  and  Y.  H.  Pao,  "On  the 
determination  of  phase  and  group  velocities  of  dispersive  waves  in  solids," 

J.  Appl.  Phys.,  48,  pp.  4320-4327,  1978. 

The  18  specimens  and  their  thicknesses  are  listed  in  Table  1.  The  phase 
velocities  as  calculated  from  the  phase  spectrum  of  pulse  propagation  for  each 


specimen  are  shown  in  Figure  3,  and  their  values  listed  in  Table  2  for  5 
frequencies. 

The  data  for  attenuations  of  18  specimens  are  shown  in  Figure  4  for 
both  P-waves  and  S-waves.  The  data  for  P-waves  are  smoother  than  those  for 
S-waves.  Even  so,  the  data  of  the  low  end  and  high  end  of  the  spectrum 
should  be  discarded. 

The  signal  of  specimen  #4  is  not  accurate  enough  for  phase  velocity 
measurement  of  the  shear  wave.  For  shear  attenuation  measurement,  due  to 
high  attenuation  or  the  mismatching  of  impedances  of  a  few  specimens  to  the 
transducer,  we  were  unable  to  obtain  the  second  echos,  and  thus  the  attenua¬ 
tion,  for  specimens  #2,  4,  11,  13  and  17. 

For  the  frequency  range  considered,  the  attenuation  constants  can  be 
represented  by  the  polynomial  approximation 

“  =  a0  +  a]f  +  a2f  2 

The  coefficients  aQ,  a-|  and  ^  are  listed  in  Table  3  for  all  materials. 


Table  1.  Dimension  of  18  Polymeric  Materials 


No.  Specimen  Thickness  (cm) 


POLYSULFONE 

1.28 

UV  curable.  Acrylate  based,  polymer 

0.694 

White  Delryn 

1.28 

Polyurethane  (Conap) 

0.715 

High  density  Polyethylene 

0.947 

Noryl 

1.51 

FEP  100 

0.803 

TEFZEL 

1.303 

REXOLITE 

1.269 

POLYSTYRENE 

1.273 

POLYPROPYLENE  7823 

0.65 

POLYYINYL  -  CHLORIDE 

1.26 

HYTREL 

0.644 

NYLON  101 

0.68 

ACRYLIC 

1.286 

TEE  (8.  6.) 

0.949 

DELRYN  157 

1.257 

POLYCARBONATE 


1.26 


Frequency 


Specimen 
#1  ' 

2  MHz 

osn 

4  MHz 
0178 

6  MHz 
05S3 

8  MHz 
0?S6 

10  MHz 

05 

2 

0.2233 

0.2265 

0.2284 

0.2296 

0.2300 

3 

0.2442 

0.2460 

0.2474 

0.2488 

0.2497 

4 

0.1675 

0.1696 

0.1708 

0.1715 

0.1658 

5 

0.2391 

0.2400 

0.2406 

0.2409 

0.2409 

6 

0.2694 

0.2701 

0.2707 

0.2709 

0.2709 

7 

0.1271 

0.1297 

0.1291 

0.1261 

0.1259 

8 

0.1647 

0.1671 

0.1660 

0.1658 

0.1647 

9 

0.2348 

0.2351 

0.2353 

0.2354 

0.2352 

10 

0.2314 

0.2316 

0.2319 

0.2319 

0.2318 

11 

0.2498 

0.2544 

0.2569 

0.2583 

0.2591 

12 

0.2343 

0.2353 

0.2360 

0.2363 

0.2364 

13 

0.2148 

0.2180 

0.2199 

0.2208 

0.2214 

14 

0.2633 

0.2654 

0.2668 

0.2675 

0.2676 

15 

0.2749 

0.2754 

0.2759 

0.2761 

0.2760 

16 

0.1378 

0.1393 

0.1394 

0.1377 

0.1378 

17 

0.2440 

0.2436 

0.2447 

0.2457 

0.2461 

18 

0.2303 

0.2306 

0.2309 

0.2310 

0.2310 

18 


6  MHz 


8  MHz 


Specimen 

2  MHz 

Frequency 

4  MHz 

6  MHz 

8  MHz 

1 

0.09303 

0.09340 

0.09366 

0.09461 

2 

0.08879 

0.08529 

0.08351 

0.08542 

3 

0.1036 

0.1053 

0.1033 

0.1034 

4 

5 

0.09240 

0.09406 

0.09322 

0.09494 

6 

0.1212 

0.1216 

0.1219 

0.1220 

7 

0.03676 

0.03459 

0.03541 

0.03541 

8 

0.05819 

0.05568 

0.05506 

0.05411 

9 

0.1167 

0.1166 

0.1166 

0.1184 

10 

0.1092 

0.1118 

0.1126 

0.1142 

11 

0.2056 

0.2199 

0.2071 

0.2055 

12 

0.1095 

0.1101 

0.1104 

0.1085 

13 

0.07787 

0.07550 

0.07673 

0.07485 

14 

0.1025 

0.1072 

0.1093 

0.1101 

15 

0.1402 

0.1404 

0.1407 

0.1428 

16 

0.05007 

0.04817 

0.04832 

0.04752 

17 

0.1083 

0.1092 

0.1070 

0.1053 

18 

0.1091 

0.1117 

0.1127 

0.1143 

3-7 
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Table  3.  Quadratic  Functions  for  Attenuation  Coefficients 


a  =  aQ  +  a-jf  +  agT  (np/cm) 


h  * 


(B)  S-waves 


Schematic  Diagram  of  Experimental  Set  Up 
Typical  Transient  Pulses  and  Their  Spectra 

(A)  First  Echo 

(B)  Second  Echo 

(C)  Spectra  for  Two  Echos 
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Attenuation  Coefficients 
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Velocities  -  Specimen  No. 
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Figure  3. 
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Figure  3.  Phase  Velocities  -  Specimen  No.  11 
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Velocities  -  Specimen 
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Attenuation  Coefficients  -  Specimen  No. 
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Attenuation  Coefficients  -  Specimen  No. 
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Attenuation  Coefficients  -  Specimen  No. 
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Attenuation  Coefficients  -  Specimen  No.  17 
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Attenuation  Coefficients  -  Specimen  No.  18 
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